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ABSTRACT 

We recently discussed an Eigenvector 1 (El) parameter space that provides 
optimal discrimination between the principal classes of broad-line Active Galactic 
Nuclei (AGN). In this paper we begin a search for the most important physical 
parameters that are likely to govern correlations and data point distribution in 
El space. We focus on the principal optical parameter plane involving the width 
of the H/3 broad component [FWHM(H/3bc)] an d the equivalent width ratio be- 
tween the Fe II blend at A4570 and H/?bc- We show that the observed correlation 
for radio-quiet sources can be accounted for if it is primarily driven by the ra- 
tio of AGN luminosity to black hole mass (L/M oc Eddington ratio) convolved 
with source orientation. L/M apparently drives the radio-quiet correlation only 
for FWHM(H/3) < 4000 km s" 1 which includes Narrow Line Seyfert 1 (NLSyl) 
galaxies and can be said to define an AGN "main sequence." Source orientation 
plays an increasingly important role as FWHM(H/?bc) increases. We also argue 
that AGN lying outside the radio-quiet "main sequence," and specifically those 
with optical Fe II much stronger than expected for a given FWHM(H/3bc), may 
all be Broad Absorption Line QSOs. 

Subject headings: quasars: emission lines - quasars: general - line: formation - 
line: profiles 
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1. Introduction 

We have recently identified a correlation space for broad line AGN involving: (1) Balmer 
line width (FWHM(H/3bc)), (2) relative strength of optical Fe ll op t and H/3 emission lines (the 
equivalent width ratio of Fell opt emission in the range 4435-4685 A and the broad component 
of H/3: W(Fell A4570)/W(H/3 B c)=RfcIi), (3) soft X-ray photon index r soft and (4) CivA1549 
broad line profile centroid shift (Sulentic et al. 2000a,b). In simplest terms the parameters 
can be said to measure: (1) the broad line region (BLR) velocity dispersion, (2) the relative 
strengths of low ionization lines that are thought to arise in the same structure, (3) the 
strength of a (thermal) soft X-ray photon excess, and (4) the amplitude of systematic radial 
motions in the high ionization gas. We call this parameter space "Eigenvector 1" (hereafter 
El) reflecting its partial origin in a principal component analysis of the low redshift (z ^ 0.5) 
part of the Palomar-Green quasar sample (Boroson & Green 1992, hereafter BG92). El 
allows us to discriminate between most AGN classes that show broad emission lines (Sulentic 
et al. 2000a,b). 

El separates the majority of radio-quiet (RQ) sources from radio-loud (RL) AGN. The 
El parameter space distribution also suggests the possible existence of two RQ classes. 
Population A with FWHM(H/3bc) ~ 4000 km s^ 1 and average El parameter values: (a) 
R Fe ii~ 0.7, (b) r soft ~ 2.8 and (c) CivA1549 centroid (blue)shift ~ -800 km s _1 . Population 
B includes all of the remaining RQ AGN with FWHM(H/3 BC ) ~ 4000 km s _1 and shows 
average El parameter values: (a) RfcII ~ 0.4, (b) r so f t ~ 2.3 and, (c) CivA1549 centroid 
shift ~ km s _1 . See Sulentic et al. (2000b) for sample variance and other details. RL 
and RQ population B sources occupy a similar parameter domain in El and show a large 
number of other observational similarities (Sulentic et al. 2000c). 

After reviewing the occupancy of the optical El parameter plane (FWHM(H/9bc) vs - 
Rfcii) and defining a "main sequence" for RQ Population A, we show that the mean ionization 
level of the broad lines decrease as one goes from RQ population B to population A sources 
(§2.2). In §4, we show how this result and the occupancy of the parameter plane can be 
explained in terms of different values of the Eddington ratio (analyzed in §3.2) convolved 
with the effect of source orientation (discussed in §3.1 and 3.2). We also consider potentially 
important outlier sources in §2.1 and in 6.1. 
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2. Relevant Trends 

2.1. Definition of "Main Sequence" and "Outliers" 

Population A sources show a clear and significant correlation among the El parameters 
while RQ population B and RL sources show a larger scatter with no obvious correlation. 
RQ population B sources occupy the same El parameter domain as, especially, flat spec- 
trum (core dominated) RL sources (see Table 2 in Sulentic et al. (2000b)). It is however 
important to recognize that uncertainties are larger for most measures of RL and RQ pop- 
ulation B sources. Limited S/N and line blending limit the accuracy of Fell opt equivalent 
width measurements for sources with W(Fell A4570) ;$ 20 A (RfcII ~ 0.2). That is why 
we cannot rule out the possibility that the RQ population A correlation extends into the 
RQ population B domain. At any rate, the RL sources are preferentially found in the same 
El domain as RQ population B, and are rarely found in the domain of population A. The 
concepts of a RQ population A - population B difference and of a RQ population B - RL 
similarity were motivated by the optical parameters but they are reinforced by X-ray r so f t 
and UV CivA1549 line shift measures. These differences/similarities are true irrespective of 
the reality of a parameter space break between population A and B (Sulentic et al. 2000b). 

Figure 1 presents a schematic view of source occupancy in the FWHM(H/3bc) vs. Rp e ii 
parameter plane. We show correlation trends for: (1) our sample (Sulentic et al. 2000b), 
(2) two radio-loud samples (Brotherton 1996; Corbin 1997), and (3) a soft X-ray selected 
sample (Grupe et al. 1999). The solid lines representing (1) and (3) connect average values 
for sources in FWHM(H/? BC ) ranges 0-2000, 2000-4000 and > 4000 km s" 1 respectively. 
All samples show the same general trends with broad line, Fell op t-weak RQ population B /RL 
sources displaced towards the upper left and narrow line, Feii op t strong, population A RQ 
towards the lower right. These lines indicate a tendency for AGN to lie along a "main 
sequence" (MS). In addition to some prototype sources, we show data points for Fell opt 
strong (FIR bright) quasars (Lipari et al. 1993). The Lipari et al. (1993) sources are shown 
as a single point representing their average FWHM(H/3bc) and RfcII, with the exception of 
two objects, IRAS 0759+651 and Mark 231 which are reported individually in Fig. 1. We 
also identify all known broad absorption line (BAL) QSOs with absorption W(CivA1549) 
£ — 9 A from the Sulentic et al. (2000b) and Lipari et al. (1993) samples. For PG quasars, 
data on absorption W(CivA1549) were obtained from Brandt, Laor, & Wills (2000) except 
for PG 1351+234, for which we measured the CivA1549 absorption on the IUE spectrum 
SWP54205. 

We see that I Zw 1 (and other NLSyl), as well as several (low z) BAL QSOs, are 
located toward the high R Fe ii end of the MS. One could infer from Figure 1 that the MS may 
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extend from the broadest double-peaked RL sources to the narrowest and Fell op t strongest 
sources (like PHL 1092). It is not yet clear if these extrema should be considered outliers or 
extensions of the correlation found for the bulk of RQ Population A sources studied so far 
(see also §6). We can identify one clear class of outliers: sources occupying the upper right 
quadrant of the El plane (Fig. 1, which turn out to be all BAL QSOs. These are sources 
with FWHM(H/9bc) ~ 3000 km s _1 and RfcII ~ 1)- In general, sources with Rf c ii>> 
[FWHM(H/3bc)/(2500 km s" 1 )]" 1 - 4 (which approximately defines the upper boundary to the 
MS of Sulentic et al. (2000a)) are sources with Fell opt in excess of the mean value expected 
for a given FWHM(H/3bc)- They will be briefly considered in §6.1. 



2.2. An Ionization Decrease from population B to population A 

Table 1 shows relevant mean parameter values for RQ population A and B sources as 
well as for the extreme population A NLSyl sources. All averages and sample standard 
deviations are from Sulentic et al. (2000b) except for I(Si III]A1892)/I(Cm]A1909) which 
comes from HST data (Laor et al. 1994, 1995; Wills et al. 1999). The decrease in equivalent 
width of CivA1549 (our representative high ionization line) along with an increase in W(Fell 
A4570) (low ionization emission) suggest a systematic decrease in ionization level from RQ 
population B to population A. We interpret the data in Table 1 by considering the behavior 
of I(Si III]A1892)/I(Clll]A1909), W(CivA1549) and W(H/3 B c) as a function of the ionization 
parameter (U) and electron density n e . We compare the observed values with a grid of 
CLOUDY computations for AGN broad line emission (Korista et al. 1997). The models 
assumed a total column density N c ~ 10 23 cm -2 and a standard AGN continuum (model 3- 
23 of Korista et al. (1997)). Figure 2 shows that all trends passing from population B towards 
population A consistently suggest a decrease in U and an increase in n e (from log(U) ~ — 1 
to —1.5, log(n e ) ~ 9.5, canonical value from AGN photoionization models to log(U) ~ —2 to 
-2.5, log(n e ) ~ 11.5). Notably, the density sensitive line ratio I(Si III]A1892)/I(Clll] A1909) 
(almost independent of U for log(U) ^ —0.5) indicates that log(n c ) ~ 10.5 - 11 towards 
the NLSyl domain. 

These considerations quantify a general trend which is very appreciable in the spectra 
of AGN with different Balmer line widths (see Fig. 2 of Sulentic et al. (2000a)) and which 
has been systematically ignored in photoionization computations. It is probably the origin 
of our inability to explain some line ratios in quasar spectra (Sulentic et al. 2000b) and is 
obviously a zeroth-order result. Our earlier comparison of CivA1549bc and H(3bc properties 
motivated us to suggest that high and low ionization lines are not emitted in the same region, 
at least in population A sources (Marziani et al. 1996; Dultzin-Hacyan et al. 2000). 
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3. The Main Physical Parameters 

3.1. The Role of Orientation in El 

Why is orientation important and why, at the same time, can it not account for all of the 
phenomenology? Important evidence in favor of orientation effects in RQ AGN involves H/3 B c 
- Civ A 1549 B c profile comparisons (Marziani et al. 1996; Sulentic et al. 2000b). Some NLSyl 
sources such as I Zw 1 show a CivA1549bc profile that is almost completely blueshifted 
relative to a very narrow H/3bo This robust observational result is easily explained in terms 
of: (1) a high-ionization wind emitting CivA1549bc an d (2) an optically thick disk emitting 
H/Sbc- The disk will obscure the opposite (receding) side of the high-ionization outflow 
(Marziani et al. 1996). Not all NLSyl show large amplitude CivA1549 blueshifts but the 
currently observed range is from to —5000 km s -1 . NLSyl also show low W(CivA1549) 
(Rodriguez-Pascual et al. 1997). The rarity of large amplitude CivA1549 blueshifts points 
toward a role for orientation since CivA1549 shifts are expected to be strongly orientation 
dependent in a disk + wind scenario. 

In principle, it is possible to ascribe low W(H/3 B c), W(CivA1549) and WQOiii] A5007) 
to an orientation dependent "blue bump" that amplifies the UV/optical continua and whose 
contribution to the continua increases with decreasing inclination (Marziani et al. 1996). 
However, this assumption is rather ad hoc and can not easily account for the increase of RFeii 
in sources with narrower FWHM(H/3bc)- 111 that scenario RFeii should remain constant 
because H/3bc an d Fell opt would be similarly affected by the amplified continuum unless 
H/9bc an d Fell opt are both, but differently, anisotropic (Marziani et al. 1996). We remark 
that the observed RFeii increase appears to be mainly due to a decrease in W(H/3bc) towards 
the NLSyl domain. This has been interpreted as the effect of collisional suppression of H/5 
(Gaskell 1985) implying an increase in electron density. The correlation between n e -sensitive 
I(Si III]A1892)/I(Clll]A1909) ratio and FWHM(H/5 BC ) also cannot be explained in terms of 
orientation. Physical conditions must change significantly along the AGN main sequence. 

Further difficulties for an orientation-only hypothesis may involve the range of [Olll] A4959,5007 
luminosity observed among PG quasars (BG 92). We suggest that this issue needs reconsid- 
eration also because the [Olll] A4959,5007 emitting region in several Seyfert galaxies shows 
a bipolar structure (Falcke, Wilson, & Simpson 1998), a result that suggests a strong orien- 
tation dependence. 
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3.2. The Eddington Ratio: A Parameter Affecting BLR Physical Conditions 

El shows us that the diversity of AGN properties can be organized on the basis of a 
set of parameters involving: r so ft, FWHM(H/?bc) ; RfcII and CivA1549 line shift. Apart 
from the reasons outlined above, it is reasonable to infer that orientation alone would be 
unlikely to account for the El correlations because we would expect to find an overlapping 
domain for all RQ and RL AGN. A distinct RL sequence, if it exists, is apparently displaced 
from the RQ one. El suggests that it begins near FWHM(H / 9 B c)~4000 km s" 1 (Sulentic 
et al. 2000a,b). Our BG92 dominated sample suggests that RQ sources become rare above 
FWHM(H/?bc) ~6000 km s^ 1 while RL are common from 4000 to, at least, 8000 km s _1 . 
Support for a RQ-RL displacement comes from the detection of a strong Ha line in BLLAC 
with FWHM(Ha)~4000 km s _1 (Corbett et al. 2000) which should be a near pole-on RL 
source. The presence of the soft X-ray excess as one of the principal correlates in El suggests 
that the Eddington ratio (i.e., the ratio L/M) may be the most important physical parameter 
driving El, as well as the main factor accounting for the RQ and RL sequence displacement 
(BG 92; Pounds et al. 1995; Boiler et al. 1996; Sulentic et al. 2000a). 



3.3. A Correlation Between FWHM(Hftc) and L/M? 

We often refer to the L/M ratio rather than the Eddington ratio (ratio between the 
bolometric and Eddington luminosities = dimensionless accretion rate hi oc L/M) because 
we rely, in this context, on independent observational measurements for both L and M. 
X-ray variability determination of M for a few AGN (Czerny et al. 2001) suggests an anti- 
correlation between L/M and FWHM(H/9bc)- Figure 3 shows the best current observational 
evidence for the anti-correlation between dimensionless accretion rate and FWHM(H/9bc)- It 
involves the sources with most accurate reverberation and/or X-ray variability based black 
hole mass determination. Sources like NGC 4051 and 4151 show a transient broad line 
component, which at some epochs is completely absent in unpolarized light (see e.g., Ulrich, 
Maraschi & Urry (1997) and references therein). FWHM(H/5bc) and X-ray properties for 
such sources may not allow a reliable mass estimate that is comparable with AGN showing 
a more typical variability pattern. At the other extreme, the two highest m radiators, 
PHL 1092, and IRAS 13224—3809 show extreme optical properties (see §6 for a possible 
interpretation). The Czerny et al. (2001) dataset is obviously biased towards objects that 
show large amplitude X-ray variability. Only six sources remain in the X-ray sample if we 
omit them which is too few for a reliable determination of the correlation coefficient. The 
results shown in Figure 3 are suggestive with the X-ray based points showing a correlation at 
a 2cr confidence level. If only the six sources are considered, Pearson's correlation coefficient 
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ip is m —0.85 with a probability P« 0.04 that uncorrelated points could give rise to the 
computed r P . If all data points are taken into account, r P pa —0.75 (P« 0.02). 

Other approaches favor a strong correlation between m and FWHM(H/3 BC ). Modeling of 
a radiation pressure-driven wind (Nicastro 2000; Witt et al. 1997) predicts the relationship 
between m and FWHM(H/3 B c) (shown as the thick line of Fig. 3). A rather strong correlation 
between L/M and FWHM(H/3 B c) also emerges using masses derived from reverberation 
mapping (Kaspi et al. (2000)) (filled circles and squares in Fig. 3). If all RQ and RL 
objects are considered except the outliers NGC 3227 (another transient H/5bc object) and 
NGC 4051, the correlation coefficient is rp pa -0.75 (P~ 2xl0~ 4 ). We obtain the following 
functional relationship by a robust fitting technique (e.g., Press, Teukolsky, Vetterling, & 
Flannery (1992)): 

(^)o » 6.2 x 10 3 FWHM l000 (H(3 BC )- 2 , (1) 

where L A is the specific luminosity at ~ 5000 A and FWHM(H/3 BC ) is expressed in units 
of 1000 km s _1 . In order to transform to m we assume a constant bolometric correction 
~ —2.5, appropriate for the typical AGN continuum as parameterized by Mathews and 
Ferland (1987). The m-FWHM(H/3 BC ) best fit is shown in Fig, 3 as a thin solid line. 

A FWHM(H/9bc) ~~ L/M correlation in the case of reverberation masses is not surprising 
since it results in part from circular arguments (i.e., FWHM(H/3 B c) is used to compute M 
from reverberation mapping data). Monte Carlo simulations were carried out assuming: (1) 
that distance r, FWHM(H/3 BC ) and luminosity are randomly distributed and uncorrelated 
in the observed ranges and (2) that M is related to the velocity dispersion by the virial 
relationship: Mocrv 2 . Observational errors for logm were assumed to contribute a Gaussian 
scatter with a pa 0.15. In approximately 4000 random trials we found a probability P £ 0.05 
that a correlation coefficient as large as 0.75 would be obtained. In order to fully circumvent 
the circularity issue we also assumed that FHWM(H/3 B c) does not correlate strongly with 
either L or M separately (as noted also by Sulentic et al. (2000a)). Optical luminosity 
appears to be an orthogonal variable with respect to the El parameters (in BG92 it is part of 
their Eigenvector 2). Mass and FWHM(H / 9 BC ) have r P ^0.48. If we simulate data points for 
which: (1) the correlation between FWHM(H/3 B c) and M falls in the range 0.43 < r P < 0.53 
and (2) mass and luminosity are not correlated, we obtain a negligible probability that an 
rp pa0.75 correlation between FWHM(H/3 B c) and logm could occur randomly. 

We conclude that the circularity inherent in the mass computation cannot fully ex- 
plain the strength of the observed L/M vs. FWHM(H/? B c) correlation for reverberation 
masses. This provides much needed support for a physical relationship between L/M and 
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FWHM(H/5bc)- Given the caveats outlined above (small numbers of X-ray determined 
masses; M(FWHM) dependence in optical reverberation masses), an independent verifica- 
tion of any fn - FWHM(H/5bc) relationship would be best derived from X-ray variability 
(Braito & Marziani, in preparation). In the present study, we assume that Eq. 1 describes 
the appropriate relation. 

Differences between the disk + wind model expectation and the linear fit described by 
Eq. I are appreciable for FWHM(H/3bc) ~ 1500 km s _1 (Fig. 3). The disk + wind model 
predicts an average difference of more than one order of magnitude between NLSyl and 
other population A sources. The assumption of constant fn for all population A yields a 
lower normalized xl (~ 3.2) than the disk + wind fit (xt ~ 4.3) (see also §6 for possible 
interpretations). 

4. Why Does Ionization Level Decrease With Increasing L/M? 

The ionization parameter can be defined as 



where Q(H) is the number of hydrogen ionizing photons and r is the distance of the 
BLR from the central continuum source. U can be rewritten in terms of L/M and M, if we 
assume: 

1. Q{H) rs fLboi/ < hv >. A typical AGN continuum as parameterized by Mathews and 
Ferland (1987) yields < u >^ 9.96 x I0 15 Hz and f « 0.54. 

2. The velocity field for the low ionization line-emitting gas is mainly rotational. We 
assume that the velocity dispersion is the square root of the mean square velocity for 
a rotating annulus: 



with FWHM(fL3 BC ) = 2.35 a. 

We also consider that the ratio I(Si III]A1892)/I(Cm]A1909) is a good density diagnostic 
(almost independent of the ionization parameter) in the density range 9.5 £ log(n e ) £ 12. 
CLOUDY (Ferland 2000) photoionization computations suggest that 



U 




(2) 




(3) 



(4) 
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Since the ratio I(Si III]A1892)/I(Cin]A1909) is directly correlated with FWHM(Hftc) (Wills 
et al. 1999), adopting the FWHM(H/5 BC ) - L/M correlation yields: 

logn e w 11.1 - 1.33 log FWHM 10Q0 (H(P) « 7.72 + - log(L bol /M) . (5) 
We can use Eq. 5 to write U in terms of L^/M: 

U = 0.26(^)-; 4 1+ ^M->, (6) 

with x = 0.67, mass in units of 10 7 M Q , and the luminosity-to mass ratio in units of 10 4 
the solar value (L/M m 1.9 ergs s _1 g" 1 ). This accounts for the somewhat counterintuitive 
result that U decreases with decreasing FWHM(H/9bc) and increasing L/M. We note that 
for M 0i 7 = 1 we get logU « —2 for FWHM(H/9bc)~ 2000 km s" 1 as expected. However a 
shallower dependence may be possible; the L/M power is constrained within 1.1 £ 1 + x £ 
1.7 for a reasonable choice of the input parameters. For instance, if we use the relationship 
between luminosity and BLR radius t B lr oc L 07 (Kaspi et al. 2000), we would obtain 
x ~ 0.1. Values of 1 + x ^ 1 are a consequence of the assumption of a Keplerian (or virial 
or similar) velocity field (i.e., v oc l/y/r). 



5. Connecting the Observational Plane to L/M and i 

The above results allow us to relate L/M to the FWHM(H/3 B c) - R-Feii plane of El. 
Determining the relationship between Rf c ii and U is not trivial since Fell opt emission is 
poorly understood. Photoionization calculations suggest a four-fold increase in total Fell opt 
emission from logU ~ —1 to —2. We assume that Rfcii scales as total Feii op t intensity 
divided by I(H/3bc) as a function of U for an average density n e ps 10 10 cm" 3 . We use the 
calculations of Korista et al. (1997) to estimate the dependence on U. The normalization 
has been chosen following Netzer (1990) with Rpdi ~ 0.25 for XogU = — 1. We note 
that using the possibly stronger correlation between I(Si III] A1892)/I(Ciu] A1909) and R Fen 
(i.e., I(Si III]A1892)/I(Clll]A1909) « 0.1 + 0.5 R Fe ii, Wills et al. (1999)) gives a consistent 
relationship without any assumption about the relationship between U and R Fe ii- Equating 
the above relationship involving I(Si III] A1892)/I(Cin] A1909) and R Fen to Eq. 4 and using 
Eq. 5 to relate n e and L/M allows us to obtain R FeII oc 0.551ogL/M. 

If low ionization lines like H/3 are emitted in a flattened configuration, then some effect 
of viewing angle is expected. In order to take into account the effect of orientation we assume 
that: 



1. the relationships employed above are valid for an average < i 30°. We write the 
FWHM(H/3 B c) dependence on i as 



FWHM(i) = FWHM(O) + AFWHM • sini, 



(7) 



where 



AFWHM = 2[FWHM(i = 30°, 



L 



) -FWHM(O)]. 




M 



2. RfcII depends on i following a sec(i) law with a ratio RfcII that may change by a factor 
1.6 (an amplitude taken from the mean Fell op t difference between lobe-dominated and 
core-dominated RL objects). 

These assumptions allow us to reproduce the parameter space covered by sources in our 
samples by assuming logM ~ 8 in the expression of U, and FWHM(O) = 500 km s" 1 . Fig. 4 
shows a grid of theoretical values superimposed on the data points of Sulentic et al. (2000a). 
Ifl + rc < 1.7 (1+a; = 1.7 is assumed in Fig. 4), then rh somewhat larger will result for 
the same RFeii, yielding however the same qualitative behavior. For instance, for x ~ 0.4, 
log(L/M) Q 4.5 corresponds to RFeii~ 1-5. 



Our calculations do not attempt to reproduce the observed point distribution, but only 
to account for the occupancy of the parameter plane, since instrumental factors and biases 
affect the distribution. The cluster of points at Rf c ii~0.2, for example, is due to limits on 
S/N and resolution. Another source of concern involves the role of selection biases in our 
El AGN sample. RQ population A sources are favored by soft X-ray (e.g., Grupe et al. 
(1999)) and optical color-based (e.g., BG92) selection techniques while RL/RQ population 
B AGN are not. The latter sources may be seriously under- represented in the RFeii vs - 
FWHM(H/5 BC ) plane. 

On the theoretical side, a distribution of masses will blur the grid, since different masses 
would deform and displace the grid horizontally. An additional source of scatter may involve 
Fe abundance. Therefore no rigorous inference can be made about individual values of L/M 
and i from the Figure 4 grid. 

The low R F eii region (RfcII ~ 0.5) of Figure 4 suggests that orientation is responsible 
for some of the population B sources. They would fall in the population A domain if viewed 
face-on. Towards the middle of Figure 4, we see that decreasing i and increasing L/M 



6. Discussion 
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apparently have a concomitant effect (decreasing i implies decreasing FWHM(H/3bc) and 
Rpeii; increasing L/M implies the same observational trends). This is also true in the domain 
of NLSyl, which should be an L/M extremum (but not necessarily of i). The concurrent 
effects of both parameters may explain why the correlations above have been found by a 
number of workers without any contradictory result. Observational prediction of some wind 
models are sensitive to both i and L/M, and account for, at least qualitatively, the CivA1549 
shift amplitude distribution as considered by Sulentic et al. (2000b). 

Radiation pressure driven wind models (Nicastro 2000; Witt et al. 1997) predict a 
decrease of rh with increasing FWHM(H/5bc)- However, the model by Nicastro (2000) 
predicts highly super-Eddington accretion for NLSyl. Existing evidence is still sparse, but 
neither a dynamical mass determination nor X-ray mass estimate supports this prediction. 
Rather, NLSyl as the bulk of population A sources, seem constrained within 0.3^m~ 1.0 
(e.g., Puchnarewicz et al. (2001) constrain 0.3 ^ rh ^ 0.7 for RE J1034+396; a similar 
result is inferred for Akn 564 (Comastri et al. 2001). Laor (2000a) also infers rh £ 0.3 
for NLSyl from X-ray variability. RQ AGN with reverberation mapping mass estimates 
(excluding NGC 4051 and 4151) of Fig. 3 show that there is no strong discontinuity between 
NLSyl and the rest of population A (i.e., sources with 2000 km s" 1 < FWHM(H/5 BC ) < 4000 
km s -1 ). The rh values become significantly different if population A and population B are 
compared with a K-S test (the difference is also appreciable in Fig. 3). The K-S test applied 
to population A (13 sources) and population B (11 sources) yields a probability P^lxl0~ 3 
that the rh values are drawn from the same parent population. 

Most NLSyl sources radiate at logm ss 0.0, as do a sizeable fraction of the sources 
with 2000 km s" 1 < FWHM(H/3 B c) < 4000 km s _1 supporting the identification of a unique 
population up to at least FWHM(H/3 B c)~ 3500 km s _1 . PHL 1092 may be radiating at 
rh > 10 (IRAS 13224-3809 may be another case). The extreme location of PHL 1092 in 
Fig. 1 and in Fig. 3 is consistent with both an exceptionally high value of rh (>> 1) and 
a pole-on orientation (also suggested by the strong X-ray variability (Forster & Halpern 
1996)). Sources like PHL 1092 may therefore be intrinsically rare even if not peculiar in a 
strict sense. 

Recent results suggest a clear dichotomy between RL and RQ AGN in terms of black hole 
mass with RL AGN having a systematically larger black hole mass (Laor 2000b). According 
to our considerations, the optical properties of AGN should be largely transparent to black 
hole mass differences (with a dependence of U on M yielding a second order effect). However, 
the probability of having low L /M is obviously favored for large masses. This is in agreement 
with the upwardly displaced location of RL AGN in the RfcII - FWHM(H/3 BC ) plot with 
respect to radio quiet AGN. Also, in the idealized case of a sample where L constant, there 
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could be a sequence of increasing mass from NLSyl to population A and then to population 
B, with RL AGN hosting the most massive black holes. 

6.1. On the Nature of the Outliers 

A major difference we are able to identify between MS and outlying BAL QSOs is related 
to their far IR spectral index a (f u oc v~ a ) between 25 and 60 /im. The index a(25,60) is 
1.47 and 1.32 for Mrk 231, and 0759+651 respectively. Consistently, from ISO data (Haas 
et al. 2000), a(25,60)^ 1.13 for PG 0043+039 (PG 1351+236 has been detected only at 
60/im). In both Fig. 1 and in Fig. 4 all the BAL QSOs known to us from Sulentic et al. 
(2000a) and Lipari et al. (1993) with W(CivA1549) ^ — 9 A in absorption are reported. 
The three MS BAL QSOs for which there are reliable IRAS or ISO data (PG 1411+442, 
PG 1700+518, and PG 1001+054) have a(25,60)^ -0.24, 0.56, -0.33 respectively. In the 
remaining two cases (PG 1004+130 and PG 2112+059) it is not possible to compute the 
a(25,60), but detection at 12/zm and no detection at 60 /im argues against large a(25,60). 

The most straightforward interpretation of this difference is a continuum steeply rising 
toward the far IR due to a significant contribution from circum-nuclear star formation in the 
outlying BAL QSOs (in line with the analysis of Haas et al. (2000) of the far IR spectral 
shapes of PG quasars). This contribution may not be dominant in the MS BAL QSOs. 
For both 0759+651 and Mrk 231, several lines of evidence suggest the presence of a strong 
circumnuclear starburst affecting the integrated broad line spectrum of these AGN (Taylor 
et al. 1999; Lipari 1994). If we consider the general population of AGN in a diagram 
rest-frame W(Fell A4570) vs. W(H/9bc), we see that Fell opt strong quasars (their average is 
shown in Fig. 1, Lipari et al. (1993)) define a boundary with W(Fell A4570)^ W(H / 9 B c), 
while the general population of AGN fills the area with W(Fell A4570) < W(H/3 B c)- Mrk 
231 and 0759+651 remain outliers: they show W(Fell A4570) significantly larger than that 
expected from W(R(3 BC ). This may indicate a significant overproduction of Fell opt due to an 
additional excitation mechanism (i.e., shocks), possibly associated with strong circumnuclear 
star formation. 



7. Conclusion 

Our work (Sulentic et al. (2000a) and references therein) shows that optical Fell op t 
emission is a fundamental parameter in AGN correlation studies. This result underlies the 
need for more sophisticated models for the production of Fell opt (e.g., Verner et al. (1999)). 
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It also points out the need for much higher S/N spectroscopic observations of many sources 
with weaker Fell opt emission in order to clarify the upper left part of El. In this paper we 
have attempted to explain, semi- quantitatively but self- consistently, the diversity among RQ 
AGN. Our attempt has been based upon the assumption that two of their most prominent 
emission features, H/3bc and Fell op t, are influenced primarily by source orientation and the 
Eddington ratio. The results have allowed us to reinforce the distinction between two RQ 
AGN populations and to tentatively identify those AGN that may be peculiar in a statistical 
and phenomenological sense including some BAL QSOs. 
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Table 1. Emission Line Mean Parameter Values for Radio Quiet AGN Populations 





W(H^bc) 

[A] 


W(Feii A4570) 

[A] 


R-Fell 


I(SiIIl]A1892)/I(Cin]A1909) 


W(CivA1549) 

[A] 


NLSyl 


72±28 


54±17 


0.8±0.3 


0.53±0.16 


41±19 


Pop A 


95±39 


60±16 


0.5±0.3 


0.43±0.19 


50±25 


Pop B 


104±34 


31±17 


0.3±0.1 


0.23±0.11 


107±85 
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Fig. 1. — The parameter plane FWHM(H/?bc) vs. RfcII- The two solid lines mark the 
average positions of the sample of Grupe et al. (1999) and Sulentic et al. (2000b). The 
dashed line traces the average loci of core-dominated (CD) and lobe-dominated (LD) RL 
AGN. Open squares identify BAL QSOs. The five not labeled objects are: PKS 1004+130, 
PG 1411+442, PG 2112+059, PG1011+054, PG1700+518 (in order of increasing R FcII ). 
Other relevant sources (see text) are also plotted. The broken scale allows for the inclusion 
of three wide-separation double-peaked RL AGN, and of PHL 1092. 
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Fig. 2.— Behavior of the ratio I(Si III] A1892)/I(Clll] A1909), of normalized W(CivA1549), 
and of normalized W(H/3) as a function of electron density n e and ionization parameter 
computed with cloudy by Korista et al. (1997). Arrow tails are placed roughly at the 
value of (n e , U) expected for Population B sources, arrow head at that for Population A. 
Equivalent width normalization is by continuum at 912 A (see Korista et al. (1997) for 
further details). 
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Fig. 3. — Relationship between dimensionless accretion rate and FWHM(H/3bc)- Filled 
circles and squares are RQ and RL AGN, respectively, with reverberation mapping mass 
estimates from Kaspi et al. (2000). Ringed filled circles label AGN with X-ray variability 
mass determination (Czerny et al. 2001). The thick line marks the prediction of a disk + 
wind model (Nicastro 2000). The thin line is the best fit employing a robust techniques 
for all data points of Kaspi et al. (2000); the thin dotted line represents the same with the 
exclusion of radio loud AGN. The vertical lines mark the boundaries of NLSyl (dotted) and 
population A nuclei (dot-dashed). 
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Fig. 4. — The FWHM(H/?bc) vs. RfcII with superimposed a grid of theoretical values as a 
function of i (10° <i< 40°) and L/M, expressed in solar values, for 3.1 < log < 4.5(=> 
logm 0), at steps of log = 0.1. A value of logM ~ 8 in solar units has been assumed 
for U. Data points are from Sulentic et al. (2000b); as for Fig. 3, filled circles represent RQ, 
and filled squares RL AGN. Data from Lipari et al. (1993) are shown as an average data 
point for objects clustering at Rf c ii~ 1-5 and FWHM(H/3 B c) ~ 2000. Open squares mark 
BAL QSOs; they are traced around a filled circle if the object was present in the sample of 
Sulentic et al. (2000b). The horizontal lines mark the boundaries of NLSyl (dotted) and 
population A nuclei (dot-dashed). The error bars in the upper right corners of the Figure 
are typical 2a confidence level errors for a data-point at FWHM(H/5 BC ) ~ 4000 km s -1 , and 
R Fe ii~ 0.5. 



